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We show that oscillations of a homogeneous medium with constant material coefficients produce 
pairs of photons. Classical analysis of an oscillating medium reveals regions of parametric resonance 
where the electromagnetic waves are exponentially amplified. The quantum counterpart of para- 
metric resonance is an exponentially growing number of photons in the same parameter regions. 
This process may be viewed as another manifestation of the dynamical Casimir effect. However, in 
contrast to the standard dynamical Casimir effect, photon production here takes place in the entire 
volume and is not due to time dependence of the boundary conditions or material constants. 



I. INTRODUCTION 



The results presented in this paper fall under the gen- 
eral category of dynamical Casimir effects. Hundreds of 
papers have been written on this subject — a review 
by Dodonov [lj published seven years ago has 333 refer- 
ences (more recent references were added in |2j). There 
is, however, a difference between all the numerous pa- 
pers on this subject and our results. The standard no- 
tion of the dynamical Casimir effect involves regions with 
boundaries, as in the original static Casimir effect - 
typical example is an oscillating mirror. In a few pa- 
pers [1, H, [f| [f| 0, H, G3 the case of changing proper- 
ties of the medium was also considered. We show that 
an infinite medium in motion without any boundaries 
and with time-independent material constants may also 
lead to photon production. This effect is rather small, 
but it represents another manifestation of the dynamical 
Casimir effect. 

The propagation of electromagnetic waves is affected 
by the motion of the medium. The simplest example of 
this effect has been predicted by Fresnel and observed 
by Fizeau in the XIX century (see a thorough review by 
Bolotovskii and Stolyarov [IH). The Maxwell equations 
in a medium moving with velocity v retain their standard 
form: 



d t D{r,t)=VxH(r,t), 
V-£>(r,t) = 0, 

d t B(r,t) = -VxE(r,t), 
V-B(r,t)=0, 



(la) 
(lb) 
(1c) 
(Id) 



but the constitutive relations between (D, B) and 



(E, H) are modified and read 

D + v x H/c 2 = e (E + v x B) , 
B v x E/c 2 = fi(H-vxD). 



(2a) 
(2b) 
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We have assumed that the medium is homogeneous 
and isotropic, characterized in its rest frame by the ma- 
terial coefficients e and /i. The constitutive relations ([2]) 
were derived by Minkowski [l2j and also by Einstein and 
Laub [l3| . Their derivations were based on Lorentz trans- 
formations and this might create an impression that they 
are valid only for a medium moving with a constant ve- 
locity. However, in the absence of dispersion these re- 
lations are valid for an arbitrary motion of the medium 
since they have a purely algebraic form (no derivatives). 
Therefore, they must be satisfied at each space-time point 
with the medium velocity taken at this point. This obser- 
vation has already been made by Minkowski who stressed 
that his relations have a purely local character — they 
hold at each time and at each point in space. In a re- 
cent book, Hehl and Obukhov [14j use geometric argu- 
ments to prove the validity of the "Minkowski consti- 
tutive relations for arbitrarily moving matter." In the 
presence of dispersion the constitutive relations do not 
have the Minkowski form ^ because they connect field 
vectors at different times. Thus, our simple analysis can- 
not be applied in the frequency regions where the mate- 
rial constants depend significantly on the wave frequency. 
Therefore, our results hold only for frequencies that lie 
in dispersion-free regions. 

In this paper we describe the influence of a harmoni- 
cally oscillating dispersion-free medium on the propaga- 
tion of electromagnetic waves. The motion of the medium 
leads to time- dependent constitutive relations which in 
turn result in parametric resonance. The existence of 
parametric resonance is a fairly general feature, but its 
influence on the propagation of electromagnetic waves 
depends on the detailed properties of the system. In our 
case, owing to the homogeneity of the medium, full anal- 
ysis is possible and there is no need to resort to various 
simplifications (replacement of the electromagnetic field 
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by a scalar field, restriction to one spatial dimension) 
that were introduced in most studies of the dynamical 
Casimir effect. 



II. MAXWELL EQUATIONS IN A MOVING 
MEDIUM 



where 



We shall seek the solution of Eqs. (JTJ) and £ 
form of the Riemann-Silberstein (RS) vector [if 



D 

V2e 



B 

i 

w 



The Maxwell equations take on the form: 



id t F(r,t) 
V-F(r,i) 



VxG(r,i), 
0, 



| in the 

M 

(3) 



(4) 
(5) 



where 



G 



E . H 

V2m \/2Y 



(6) 



The Minkowski constitutive relations ([2]), expressed in 
terms of F and G, are 



cF + inv x F = nG 



v 

i— x G, 

c 



(7) 



where n = c^/e/I is the refractive index of the medium and 
v is the local velocity of the medium. The advantage of 
using the RS vector is the reduction of two equations to 
one equation without any loss of information — we can 
always recover the original two equations by taking the 
real and imaginary parts. Equation (J7]) can be solved for 
G leading to constitutive relations in the form 



G= C - 

n 



1 



■dx (vxF + icnF) 



(8) 



The Maxwell equations 0$ and ([5]) together with con- 
stitutive relations (|8|) form a complete set of equations 
for the Riemann-Silberstein vector F. The problem at 
hand is that of electromagnetic wave propagation in a 
medium with time-dependent constitutive relations. A 
distinctive feature of the present case is that the con- 
stitutive relations cannot be separated into electric and 
magnetic parts. After the separation of Eq. © into its 
real and imaginary parts, we see that the vector E (and 
similarly H) depends on both D and B. 

From now on we will assume that the velocity of the 
medium depends only on time, and we express it in terms 
of a dimensionless vector (3(t), 

v(t)=cf3(t)- (9) 
In this case, the Maxwell equations (@| read 
id t F = -icj{t)[/3(t)-V]F 

+ (c/n)a(i)VxF - (c/n)j(t)/3(t)xV\j3(t)-F], (10) 
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7(<) n 2 -/? 2 (i)' 
a{t) = l- 1 {t)l3 2 (t). 



(11) 
(12) 



III. SOLUTIONS OF THE MAXWELL 
EQUATIONS IN A MOVING MEDIUM 

We shall construct the general solution of the Maxwell 
equations by performing the Fourier transformation with 
respect to the spatial variables. Since the coefficients in 
Eq. pop depend only on time, it is convenient to use the 
projections of the field on plane waves, 



f{k,t) 



d 3 7 



(2tt) 3 / 2 



e- lkr F(r,t). 



(13) 



The components of the Fourier transform /(fc, t) obey 
the following set of ordinary differential equations: 

= - ic ^t)[k-P(t)}f(k, t) + -a(t)kxf(k, t) 
at n 



->y(t)kxl3(t)WYf(k,t)]. 



(14) 



The first term on the right-hand side can be eliminated 
by the following change of the overall phase: 



f(k,t) = e-^f(k,t), 

<f>(t)=4 dt' 7 (t')[k-(3(t% 
Jo 



(15) 
(16) 



and we obtain 



ffl^ = C -{a{t)k x />, t)+ 7 (t)fc x (3{t) [f3(t) ■ f(k, t)}}. 

(17) 

Out of these three equations only two are independent 
since the amplitude f(k,t) must obey the divergence con- 
dition k ■ f(k,t) = 0. Therefore, we may always elimi- 
nate one equation from this set of three equations by 
taking the projections on two independent directions in 
the plane normal to k. 

In this paper, we assume that the velocity has a con- 
stant direction, specified by a unit vector m, and only 
its length changes in time. Before discussing the general 
case we shall solve the evolution equations when the wave 
vector and the direction of motion are collinear. 



A. Propagation along the direction of velocity 

In this case, the evolution equation (fl7|) reads 

df(k,t)/dt= (c/n)a(t)kxf(k,t). (18) 
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It describes a precession of the vector / (fc, t) around the 
wave vector fc with a varying angular velocity governed 
by the time-dependent velocity of the medium. One may 
check by inspection that the solution of Eq. (jT5J) is: 

f(k,t) = /(fc,0)cosV(i) + fcX/ , (M) sinV(t), (19) 



where 



where 



Mt) = — \ dt'a{t'), 
n Jo 



(20) 



and k — \k\. Since the length of the vector f(k,t) does 
not depend on time, |/(fc,i)| 2 = |/(fc,0)| 2 , parametric 
resonance does not occur. 



B. Propagation along an arbitrary direction 

Unlike the simple solution obtained above in the spe- 
cial case, the solution in the general case cannot be given 
in terms of elementary functions. To reduce Eq. (jTTJ) to 
a more amenable form, we expand the vector / (fc, t) in 
an orthonormal basis in the plane perpendicular to fc, 

f(k,t) = m(fe)/i(fe,t) + n 2 (fc)/ 2 (fc,t). (21) 

We shall use the distinguished direction m to construct 
two unit vectors, 

, fcx(fcxm) ,, . kxm . . 
"i(*0 = 77. -, n 2 {k) = - , (22) 



kk± 



where k± = \k x m\. These vectors, together with the 
unit vector k in the k direction form an orthonormal 
basis, 

ni x fc = ri2, k x n 2 = rii, n 2 x ni = fc. (23) 

The evolution equations (|17[) rewritten in terms of 
fi(k,t) and f 2 (k,t) read 



dfi(k,t) 

dt 
df 2 (k,t) 

dt 



ua(t)f 2 (k,t), (24a) 
-UK(t)h(k,t), (24b) 



where 



n(t) = 1 - cos 2 -d j(t)/3 2 (t), 



(24c) 



lu = ck/n, and $ is the angle between the wave vector k 
and the direction of the velocity m. 

In Sec. |V] we will need a different decomposition of 
/(fc, t) that involves a complex polarization vector, 

f(k,t) = e(fc)/+(fc,t) + e*(k)f-(k,t), (25) 



e(fc) = 

/+(*,*) = 
f-(k,t) = 



ani — *c n 2 

<T- 1 / 1 (fc,t)+ i( T/ 2 (fc,t) 

V2 

a^hik^-iahjk^) 
V2 



(26a) 
(26b) 
(26c) 



At this moment we do not specify the value of the real 
parameter a. The functions f±{k,t) satisfy the following 
evolution equations: 



df+(k,t) 
dt 

df-(k,t) 
dt 

where 



= -<« {v+ (t)f+ (fc, t) - V- (t)f- (fc, *)] , (27a) 
= iw [»,+ (*)/_(*,«)-»,_(«)/+(*,«)], (27b) 



(28) 



To obtain concrete results we must specify the time 
dependence of the velocity. 



C. Motion with constant velocity 

Before analyzing accelerated motions, we shall solve 
the case of steady motion. This result will be later used as 
a reference point in the definition of the photon creation 
and annihilation operators when we will need the decom- 
position of the field operators into positive and negative 
frequency parts. 

For a constant velocity, the functions a(t) and n(t) 
become time independent. If a = (a/n) 1 ' 4 , then 7/_ 
vanishes, and Eqs. (|27| decouple. Their solutions are 



/ + (fc,<) = /+(fc,0)e— l4 , 
/_(fc,<) = /_(fc,0)e^ 4 , 



(29a) 
(29b) 



where uj\ — uj^Jqlk. 

For steady motion, the phase <f>(t) is a linear function 
of time, 



<p(t) = Wot, UJq = ck^ft cos-d. 
Thus, the complete solution of Eq. (jT^J) is 



/(M) 



f(k,t), 



f(k,t) =e(fc)/ + (fc,0)e 



e*(fc)/_(fc,0)e 



(30) 



(31) 
(32) 



This formula does not always give a decomposition into 
positive and negative frequency parts because both u)\ + 
<x>o and wi — luq may change sign. However, this can hap- 
pen only when the velocity v exceeds the speed of light 
c/n in the medium, because then luo may become larger 
than uj\. We shall restrict ourselves to moderate veloci- 
ties when this interesting phenomenon does not occur. 
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D. Harmonic oscillations 

We shall assume harmonic oscillations of the medium 
with the amplitude A and the angular frequency CI, 



(3{t) = fccos(nt), 



(33) 



where b = ACl/c. The formulas for j(t) and a(t) take on 
the forms 



7(*) 



1 



n 2 — b 2 cos 2 (Clt) ' 
a(t) = 1 - j {t)b 2 cos 2 (fit). 

For harmonic oscillations, the phase (jTHJ) is 



ck cos d 



1 



Vn 2 - 6 2 



arctan 



sin(Qi) 



(34) 
(35) 

(36) 



It vanishes when the wave vector k is perpendicular to 
the direction of oscillations. 

Equations ([24)) can, of course, be solved numerically. 
However, to get a better understanding of the general 
properties of the solutions, we shall reduce the two equa- 
tions to a single equation for a parametrically driven me- 
chanical harmonic oscillator. This will lead us, in the 
small-6 approximation, to the Mathieu equation whose 
properties can easily be explored with the use of MATH- 
EMATICA. The appearance of the Mathieu equation in 
the dynamical Casimir effect has been noticed before 
[l7l [lH in the case of an oscillating boundary. 



E. Mathieu equation 

To obtain the Mathieu equation, we make the substi- 
tutions 



/l(M)=«V«(*)Q(M) 

1 d 



/ 2 (M) 



Vaft)Q(M) 



a(t) dt 

and we obtain the following equation for Q(k,t): 
d 2 Q(k,t) 



dt 2 



h(k,t)Q(k,t) = 0, 



where 



h(k,t) = u) 2 a(t)n(t) - a/oW 



dt 2 \y^{t) 



(37a) 
(37b) 

(38) 
(39) 



After the decomposition of h(k, t) into the Fourier se- 
ries, E q. (13"51) takes on the form of the Hill equation 



d 2 Q(r) 
dr 2 



6o + 2^T9i cos(2/r) ] Q(r) = 0, (40) 



where r = Clt. We have omitted the argument k in Q(t) 
but this dependence is still there through the parame- 
ters 9i . Under normal conditions the parameter b is very 
small because the velocity of the medium is much smaller 
than c. We are, therefore, fully justified (see Fig. [2] for 
a graphic demonstration) to resort to the nonrelativistic 
approximation and keep only the hrst correction in the 
expansion of h(k,t) with respect to b 2 . In this approxi- 
mation Eq. (|38|) becomes the Mathieu equation, 



+ [a - 2 g cos(2r)] Q(t) = 0, (41) 



d 2 Q(r) 
dr 2 



whose properties are known in great detail [22], |23|, |2J] . 
In our case the Mathieu parameters a and q have the 
values 



1 



b 2 n 2 



2 n 2 



(l + cos 2 i9) 



to \ 2 b 2 n 2 — 1 , „, b 2 n 2 — 1 

-1 ^(l + cos 2 tf) 



(42a) 



(42b) 



.CI J 4 n 2 v ' 2 n 2 

The most important property of the Mathieu equation 
is the appearance of alternating stable and unstable re- 
gions — a clear signature of parametric resonance. In 
the unstable regions the solutions grow exponentially. 



IV. PARAMETRIC RESONANCE 

The unstable regions in parametric resonance lie 
around the values of u equal to multiples of CI. The 
location of these regions for the Mathieu equation is well 
known (see, for example, [12]). The two lowest unstable 
regions in the plane of the parameters a and q are shown 
in Fig. [TJ A typical behavior of the amplitude in the 
region of parametric resonance is shown in Fig. [21 The 
two almost coinciding curves are the best proof that the 
approximation of the solutions of Eqs. ([2~4")l by the Math- 
ieu functions is very good even for relatively large b. To 
describe the properties of the electromagnetic field in a 
given medium we shall vary the wave parameters lo and 

keeping n and b fixed. Both a and q depend linearly 
on J 2 . Therefore, all the points in the (a,q) plane that 
differ only in the values of ui lie on a straight line defined 
by the equation 



b 2 n 2 - 1 



■(1 + cos 2 tf)a 



= 1 



b 2 n 2 - 1 



1 - 



b 2 n 2 - 1 



(l + cos 2 tf)^ , (43) 



as shown in Fig. [T] Thus, parametric resonance occurs 
for those values of a and q where the dashed line in this 
figure crosses the shaded areas. With increasing u> the 
line crosses the consecutive unstable regions. 

According to the Floquet theory (see, for example, [13, 
HE Hll ) every solution of the Hill equation ([4"0)) has the 
form 



z=i 



Q(r) 



e WT P{r) + e 



~P(-t), 



(44) 
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FIG. 1: Unstable regions (shaded areas) in the (a, q) plane. 
The dashed line, defined by formula (|43p . corresponds to the 
values n — 2, b — 0.3, and t? = tt/2. The same values are 
chosen in Figs. [2}{4] and [6] 



where P(t) is a periodic function with the period equal 
to 7r. The essential information about parametric reso- 
nance is contained in the characteristic exponent v. In 
the region of stable oscillations the characteristic expo- 
nent v is real. The imaginary part of v in an unstable 
region determines the rate of the exponential growth of 
solutions with time. 

All numerical results presented in this paper are based 
directly on the solutions of Eqs. ([M]) and [[27]), Their 
connection with the characteristic exponent appearing 
in the Floquet representation (I44| is explained in the 
Appendix. 

It is the general property of all Hill equations with 
real coefficients (see, for example, [2(|) that in each un- 
stable region the real part of v is constant, equal to a 
natural number. Therefore, the plot of the real part 
of v will exhibit characteristic steps, as seen in Figs. [3] 
and 31 Unfortunately, the consecutive regions of para- 
metric resonance become narrower and narrower. Each 
time we move to the next higher resonance region the 
step shrinks by roughly a factor of b 2 and the imaginary 
part of v diminishes at the same rate. Thus, the low- 
est resonance plays the dominant role. For the lowest 
resonance, the values of the imaginary part of the char- 
acteristic exponent with varying b are depicted in Fig. [5] 
Their maximal values can be approximated by the for- 
mula Imi/Ri 0.076 2 + 0.136 4 . 




FIG. 2: In this plot we see two curves that lie almost exactly 
on top of each other. One represents the numerical solution of 
Eqs. (|24p and the other the analytic solution constructed from 
the corresponding Mathieu function for the same parameter 
values as in Fig. [T] We have chosen the value lo/Q, = 1.016 
for which the imaginary part of the characteristic exponent is 
largest (see Fig. [4}. The exponential growth of the amplitude 
in time is clearly seen. 




(o/a 



FIG. 3: Real part of the Mathieu exponent plotted as a 
function of ui/Q showing a tiny step at the first resonance 
(Rez/ = 1). The step at the second resonance (Rez/ = 2) is 
too small to be seen on this scale (see Fig. |4}. 



the Fourier transform of the RS vector, 

f{k, t) = e -*<*> (ni(k) uv^I Q(k, t) 



n 2 (k) 



1 d 
a(t)di 



(45) 



From the solutions Q(t) of Eq. 



we may construct 



where, for clarity, we have reintroduced the dependence 
of Q(t) on fc. The essential features are determined by 
the function Q(k,t), since the phase factor exp[— i<j)if)\ 
does not influence the discussion of parametric resonance. 

The presence of parametric resonance for electromag- 
netic waves propagating in an oscillating medium leads 
to the amplification of those waves whose wave vectors 
lie in the resonance regions. The largest effect is for the 
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FIG. 4: Behavior of the real and imaginary parts of the 
characteristic exponent v in two resonance regions (note the 
change of scale) . The figures to the left correspond to the first 
resonance (Re v = 1) and those to the right correspond to the 
second resonance (Re v = 2) . 



Im v 




FIG. 5: Plots of the imaginary part of the characteristic ex- 
ponent v in the first resonance region for consecutive values 
of 6 from 0.4 to 0.9. 



waves with the frequencies close to ft propagating in the 
direction perpendicular to the direction of oscillations. 
However, it would be difficult to achieve a significant 
amplification because it is governed by the exponential 
factor exp(0.1b 2 flt). The essential parameter b in this 
expression — the maximal velocity of the medium di- 
vided by the speed of light — is very small for realistic 
situations. It would take many periods of oscillations to 
achieve an observable amplification or significant photon 
production studied in the next section. 



V. PHOTON PRODUCTION 

The dynamical Casimir effect may be viewed as an am- 
plification of the vacuum fluctuations of the electromag- 
netic field that leads to the production of real photons. 
These effects are enhanced in the region of parametric 
resonance where the amplification grows exponentially 



with time. 

The starting point of our analysis is the replacement of 
the classical field F(r,t) by the field operator F(r,t) in 
the Heisenberg picture. Since this operator satisfies the 
same field equations as its classical counterpart, we can 
use the solutions of the field equations obtained in previ- 
ous sections. The replacement of the classical functions 
by operators leads to the following representation of the 
field operator F(r, t) 



F(r,t) 



d 3 k 



. e ik-r e -i<t>(t) 



(2tt) 3 / 2 

e(fe)/+(M)+e*(fc)/-(M) 



(46) 



where the polarization vector e(fc) is given by Eq. (|26p . 
From the equal-time canonical commutation relations 

--e tjk d k 6^(r-r'), (47) 



Fi(r,t),F](r',t) 
we obtain 

~f + (k,t),fl(k',t)] = hu6&>(k 
/_(fc,<),/!(fc',<) 
f-(k,t)Jl(k',t) 



= = 



k'), 

hwS (3 \k- k'), 
/_(*,«),/+(*',*) 



(48a) 
(48b) 
(48c) 



A. Creation and annihilation operators 



The opposite signs in the commutation relations (I48aj) 
and (|48bp may suggest that f + (k,t) plays the role of an 
annihilation operator of one type of photon while /_ (fe, t) 
plays the role of a creation operator of another type of 
photon. However, such an interpretation requires sepa- 
ration into positive and negative frequency parts. This 
is possible only for a time-independent velocity. We can 
make use of this possibility if we assume that the velocity 
is constant at the beginning and at the end of the oscilla- 
tory motion, as shown in Fig. [5^. In the regions of steady 
motion we can split the operator f(k, t) into positive and 
negative frequency parts. This enables us to define the 
incoming and outgoing creation and annihilation opera- 
tors for both types of photon. We use again the value 
of <7 = (oj/k) 1 / 4 corresponding to the initial (and final) 
velocity, 



t<0 f+(k,t) = Vhwa in (k)e- tuJlt , 

/_(*,<) = Vn^bl n (-k)e iwit (49a) 
t > t n f+(k,t) = %/^aout(fe)e-^ l( *- tN) , 

/_(fc,i) = ^6t ut (-fe)e^ l(t -* N) , (49b) 

where tjy = 2ttN/H specifies the number of periods of 
the oscillatory motion. The reversal of the sign of k in 
the second equation follows from the interpretation of 
/_ as the creation operator. The factor y/huj has been 
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extracted to secure the standard normalization of cre- 
ation and annihilation operators. In the limiting case of 
a motionless medium we obtain the standard annihilation 
operators of photons with right and left hclicity, 



Oi n (fe) = a out (fe) = a R (k), 
Sin(fe) = S ut(fe) = a L (k). 



(50a) 
(50b) 



The connection between the incoming and outgoing 
creation and annihilation operators can be obtained by 
solving the differential equations (f27|) for the operators 
f±(k,t). The time evolution from t = to t = tjy leads 
to 



flout (fe) 



fi(k, t N )a ia (k) + fi(k, t N )bl(-k) 



fl(k,t N )a in (k) + f 2 _(k,t N )bi(-k) 



(51a) 



(51b) 



where the solutions 
{fl(k,t),f 2 _(k,t)} of Eqs. 
initial conditions: 



(f2"T)) obey the following 



/|(fe,0) = l, fL{k,Q) = 0, (52a) 
fl(k,0)=0, /!(fe,0) = l. (52b) 

The mixing of the creation and annihilation operators 
(Bogoliubov transformation) is a clear signature of a two- 
mode squeezed vacuum state, accompanied by the cre- 
ation of photon pairs. In order to calculate the photon 
production rate we assume that at t = there are no 
photons. It means that the state of system \i[>) is char- 
acterized by Oin(fc)IV') = an d bi n (k)\ip) = 0. 



B. Number of photons 

In infinite space the average number of produced pho- 
tons is, of course, infinite. This infinity appears in the 
expression for the average number of outgoing photons, 
evaluated in the state \ip), as the Dirac 5 function with 
its argument equal to zero, 

MaL(*Ofiout(fcM - \fHk,t N )\ 2 S^(0), (53a) 
(#L(fe)&out(fc)IV>) = \.fl(k,t N )\ 2 S^(0). (53b) 

To make use of these formulas we have to give a physical 
meaning of the dimensionless functions \f^_(k,t)\ 2 and 
|/i(fe, t)\ 2 . First, let us note that owing to the symmetry 
of the evolution equations (|27p under the transformation 



f+(k,t) -> f*_(k,t), 
/_(fc,t) -> fl(k,t), 



(54a) 
(54b) 





FIG. 6: (a) shows the general behavior of the velocity of the 
medium (in arbitrary units) as a function of ttt. The two 
insets show magnified views of areas where the constant ve- 
locity is smoothly joined to harmonic oscillations, (b) shows 
periodic oscillations of the photon number in the stable region 
(w/Q = 1.55). (c) shows an exponential growth of the photon 
number in the region of parametric resonance (w/fi = 1.016). 



the two functions in Eqs (|5"3")) are equal (both types of 
photons are produced with the same probability). Their 
common value n(k, t), 



n(k,t) = \f 1 _(k,t)\ 2 = \f 2 + (k,t)\ 2 , 



(55) 



is the number of photons of one type in the volume h 3 
of the photon phase space. This can be shown with the 
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help of the standard procedure of artificially restricting 
the volume to a finite cube of volume V = L 3 and im- 
posing periodic boundary conditions. Then, |/l(fc,i)| 2 
and \fl(k,t)\ 2 are the numbers of photons in the dis- 
crete modes specified by the wave vector k at time t. In 
the limit, when L — ► oo, we will find that n(k,t) is the 
density of photons in phase space. In other words, the 
total number of photons riR in a dispersion-free region R 
of the phase-space is given by the formula 

»*(*) = 2 Ji^ n {k,t), (56) 

where p = hk. 

Obviously, there is a big difference between the time 
dependence of n(k,t) in stable and unstable regions. In 
the stable regions we have periodic oscillations and in the 
unstable region we have an exponential growth (Fig. [6]). 

VI. CONCLUSIONS 

We have demonstrated that, in addition to the two 
well-known mechanisms (time dependent boundary con- 
ditions or time dependent material coefficients) that lead 
to photon production through the dynamical Casimir ef- 
fect, there exists a third mechanism. The main result of 
our study is the prediction of photon production in an os- 
cillating medium without boundaries and with constant 
material coefficients. The density of photons in phase 
space given by Eq. (|55[) grows exponentially in time in 
the regions of parametric resonance. These are the same 
regions where classical electromagnetic waves are ampli- 
fied. Unlike the results given in most papers on the dy- 
namical Casimir effect, our result is exact and based on 
full Maxwell equations in three dimensions. 
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APPENDIX A: APPLICATION OF THE 
FLOQUET THEORY TO EQS. (j27) 

In order to obtain the characteristic exponents that 
identify the location of parametric resonance for a set of 
equations we need the monodromy matrix constructed 
from the fundamental solutions of these equations. Flo- 
quet theory states that the 2x2 matrix constructed from 
the fundamental solutions defined in Eqs. (f5"2"|) has the 
form 

$::;!) (A1 , 

where X(t) is a periodic function of r with period it and 
K is a constant matrix. The fundamental solutions eval- 
uated at t = 7r form the monodromy matrix M = X(n) 
which determines all essential properties of every Flo- 
quet equation (see, for example, [25|,[26j]). In particular, 
its eigenvalues fi± determine the characteristic exponent 
v. The eigenvalue equation for M is 

/i 2 -/4/j(fc,7r) + / 2 .(fc )7r )] /2 

+ [fUk,ir)fl(k,ir)-fi(k,ir)fl(k,%)] =0. (A2) 

The last term, the determinant of the monodromy ma- 
trix, is the Wronskian of Eqs. (|27|) . It is equal to 1 — its 
value at r = — because the Wronskian is a constant 
of motion. Hence, the roots of Eq. (|A2|) can be cast into 
the form 

l*±= [fl(k,n) + fl(k,n)] /2 

±%yjl- (fl(k, tt) + p_ (fe, tt)) 2 /4 = e*™, (A3) 

where v = arccos [(/.|(fe, n) + f 2 (k, tt)) /2] is, in gen- 
eral, complex. 
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